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ABSTRACT: A phenomenological analysis of K — > nn^ (Kl — > 7r + 7r~7 and 
K + — > 7r + 7r°7) with the direct emission photon is carried out beyond the leading 
order in the chiral perturbation theory. We show that the experimental evidence 
for the large photon energy dependence of the magnetic amplitude in Kl — > 7r + 7r~7 
.£h j seems to indicate an interesting consequence: vector meson dominance must be im- 

plemented at 0(p 4 ), which is not a general feature of the chiral perturbation theory. 
The phenomenology of K + — ► 7r + 7r°7 is also analyzed using the same scheme. 
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1. Introduction 

Non-leptonic kaon decays have been an important tool for studying the weak inter- 
actions [0, 0, [| £|, ^]. The radiative non-leptonic kaon decays such as the processes 
K — > 7T7T7 (i^i — » 7r + 7r _ 7 and fT + — > 7r + 7r°7) are dominated by long distance contri- 
butions, and it is not easy to match long and short distance contributions of these 
processes. In Ref. ||, an impressive improvement for the evaluation of the weak 
matrix elements in some particular processes matching long and short distance con- 
tributions is obtained in the large Nq limit. Spin-1 resonances are implemented in 
this context and important to achieve the matching. Here we somewhat complement 
their work by looking for a good low energy phenomenological description of the 
spin-1 resonances in the weak sector to study the decays of K — > tttt^. 

The total amplitude of K — > contains two kinds of contributions: the inner 
bremsstrahlung (IB) and direct emission (DE). Due to the pole in the photon energy 
the IB amplitude generally dominates unless the non-radiative one is suppressed due 
to some particular reason. This is the case of Kl — > 7r + 7r~7 and K + — > 7r + 7r°7. 
The non-radiative amplitude of the former is suppressed by CP invariance and the 
latter one is suppressed due to the AI = 1/2 rule. It is of interest to extract the DE 
amplitude of these channels in order to reveal the chiral structure of the processes. 
DE contribution can be decomposed into electric and magnetic parts in a multipole 
expansion || [?|. The available experimental evidence is consistent with a dominant 
magnetic part for the DE amplitude |§ §. So we will focus our attention on the 
magnetic part. 



In the framework of chiral perturbation theory(xPT) [p~0] , IT| , K — > nn^y has 
been analyzed previously 0, |6], 0, [L2|, [13], [14], |TJ, [L6|], and the leading order of the 
magnetic amplitudes of the processes, starting at 0(p 4 ), appear as a constant with 
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two kinds of contributions: i) the reducible type from Wess-Zumino-Witten action 
and 0(p 2 ) weak lagrangian £ 2 S=1 ; ii) the local type from 0(p A ) weak lagrangian 
C^ s=1 . However, experimental analysis has found a clear and large dependence on 
the photon energy in K L — > 7r + 7r~7 || [I7J. In the K + case, the energy dependence 
of the magnetic amplitude has not been observed yet. In order to explain the photon 
energy dependence in x?T, the theoretical analysis has to be beyond the leading 
order 0, [15], [16| . A complete 0(p 6 ) magnetic amplitude of K — > tttt^ generally may 
be useful in the future but not now since some unknown parameters have to be 
introduced, which in fact makes the prediction impossible. 

The recent direct measurement of the Kl — > 7r + 7r~7 DE form-factor by KTeV 



Collaboration [17] clearly indicates a vector meson dominance (VMD) form-factor 



T= m* + ^ 2 ' (L1) 

m K lm K 

1 T ^ 2~-^7 

rriy my ' 

where A\ and A 2 are constants with A\j A 2 = —1.243 ± 0.057 , and E* is the photon 
energy in the rest frame. Eq. (|1.1|) gives the best \ 2 f° r a single-parameter 
fit (x 2 /DOF is 38.8/27), compared with the linear slope fit (x 2 /DOF is 43.2/27), 
and two-parameter quadratic slopes fit (x 2 /DOF is 37.6/26). Therefore, this mea- 
surement seems to indicate that VMD should be implemented at 0(p 4 ) instead of 
being at 0(p 6 ) for this decay. This point is not well understood currently within 
xPT. The purpose of the present paper is to understand it using a phenomenological 
description. Also, we extend our analysis to the decay K + — > 7r + 7r°7. 

In Section 2, we remind briefly the kinematics of K — > mvy. In Section 3, we 
carry out a phenomenological analysis of K — ► 7r7r7 using chiral lagrangian plus 
VMD. The results are summarized in Section 4. 



2. Kinematics 

The general invariant amplitude of K — > tcttj can be defined as follows 0, IT3 



A[K{p) - 7r 1 (p 1 )7r 2 (p 2 ) 7 ( g ,e)] = e^q)M^q, Pl ,p 2 ), (2.1) 

where e M (g) is the photon polarization and M M is decomposed into an electric E and 
a magnetic M amplitudes as 

M» = ^-\pvqV2„ ~ Prqpi,] + ^^-e^rfjftqP , (2.2) 



with 



Q-Pi / ■ i o\ P-1 

— =-,(« = 1,2), Z 3 = — Z 3 = Z!+Z 2 . 

m z K m z K 
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The double differential rate for the unpolarized photon is 
d 2 T 



m k 



dzidz 2 (47r) 



(1^)1 + \M( Zi )\ 2 )[ Zl z 2 (l - 22 3 -r{- r 2 ) - rfz 2 2 - r'zf], (2.3) 



where r« = m^Jmx- 

In Kl — > 7r + 7r~7, the most useful variables are: (i) the photon energy in the 
kaon rest frame E* and (ii) the angle 9 between the photon and 7r + momenta in the 
di-pion rest frame. The relations between E*, 9 and the Zj are: 



E^_ 

m K ' 



z± 



E* 

7 (1=f/?cos0), 



2rriK 



(2.4) 



where j3 = J 1 — 4m^/(m 2 K — 2mxE*). Then the differential rate is 



d 2 T 



(E;fp 3 



dE*dcos9 5127r 3 m 3 K 



' x 2E*; 

m K . 



sm 2 9(\E\ 2 + \M\ 



(2.5) 



For K + — > 7r + 7T°7, three photons will be detected in the measurement, so it is 
more useful to study the differential rate as a function of: (i) the charged pion kinetic 
energy in the K + rest frame T*, and (ii) W 2 = (q-pK)(q • P+)/( m2 r+ m lr) @- These 
two variables are related to the z,- t by 



2m\ 



Z3Z+ 



mi 



m\ 



[m 2 K + ml+ - m 2 - 2m K m 7T + - 2m K T*), 



+ -W 2 . 



(2.6) 
(2.7) 



The advantage of using these variables is that, through the W 2 dependence, one 
can easily disentangle the different contributions of the IB, DE amplitudes, and 
interference term between IB and DE 

8 2 T d 2 T, 
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where A = A(K + -> 7r + 7r°) 



3. Analysis 

The leading order magnetic amplitudes of K — > ttttj start at 0(p 4 ) in xPT||, [HJ 

(4) _ eG 8 m 3 K 



M 



.1/ 



(4) 



L ~ 2n 2 F 
eG s m 3 K 

4tt 2 F 



(a 2 + 2a 4 ), 



[2 + 3(203-02)]. 



(3.1) 
(3.2) 



3 



The subscripts L and + denote Kl — > 7r + 7r~7 and K + — > 7r + 7r°7 respectively. The 
Oj's parts of the above amplitudes come from the local weak lagrangian C^ s=1 . The 
first term in is the reducible type contribution. Due to the Gell-Mann-Okubo 
mass relation, the reducible magnetic amplitude of Kl — > 7r + 7r~7 generated from 
Kl — tt°{i]8) mixing vanishes at 0(p 4 ). However, when rf is included, thus rj — rj 



mixing is considered, there is a reducible amplitude called F\ term in Refs. |4], |XG 
which is therefore at 0(p 6 ). 

The experimental analysis of Kl — > 7r + 7r~7 using the VMD form-factor parame- 
terization by KTeV Collaboration |T7| indicates that VMD must be implemented at 



0(p 4 ). This means that the couplings Oj's, £=1, 2, 3, 4 (or in terms of iVj's, 
29, 30, 31) should get the contribution from the vector resonance exchange. 

7 7T-(7r°) 7T-(7r°) 
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V \ p V 

K L (K+) 7T+ K L (K+) 7T+ 

(a) (b) 

Figure 1: Diagrams contributing to the indirect VMD magnetic amplitude of Kl —> 
7r + 7r~7 or K+ — > 7r + 7r°7. The diamond in the external legs denotes Kl — 7r°,rjs or K+ 
mixing. The black circle denotes the strong/electromagnetic vertex. The crossed diagram 
7r + <-> 7r _ (7r°) should be considered in (b). 
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-(tt°) 7r-(7r°) 7 





K L (X+) 7T+ K L (K+) 7T+ 

(a) (b) 

Figure 2: Diagrams contributing to the direct VMD magnetic amplitude of Kl — ► 7t + 7t _ 7 
or if + — > 7r + 7r°7. The empty box denotes the vertex generated by eq. ( |3.3[ ), and the 
black circle denotes the strong/electromagnetic vertex. The crossed diagram 7r + «-> 7r~(7r ) 
should be considered in (b). 

VMD can be introduced into the effective lagrangian phenomenologically, and 



4 



there are two different kinds of vector resonance exchange contributions to the non- 
leptonic radiative kaon decays |18|, [L9| : 



(i) Vector resonance exchange between strong/electromagnetic vertices with a 
weak transition in an external leg, as shown in Fig. 1, which are usually called as 
indirect transitions. The amplitude from this kind of transition is of reducible type, 
and vanishes at 0(p 4 ), which therefore contributes first at 0(p 6 ). 

(ii) The direct weak transitions are those where the weak vertices involving the 
vector resonances are present, which could contribute to the couplings a^'s. But this 
is NOT a general feature of VMD. Indeed the antisymmetric tensor realization of 
vector resonances does not generate any vector exchange contributions to K — > tttt^ 
at 0(p A ). However, this is not the case for the other realization approaches such as 
massive Yang-Mills, hidden local symmetry, and conventional vector formulations. It 
has been pointed out in Ref. |Tj| that, for the odd- intrinsic parity operator relevant in 



the V — > P7, the antisymmetric tensor formulation would give contributions starting 
at 0(p 4 ) while QCD requires an explicit 0(p 3 ) term given by the conventional vector 
formulation. As already realized in |[20|| , using the conventional vector formulation, 



there are 0(p 4 ) VMD contributions generated by the following operators through 
the direct weak transition (see Fig. 2): 

C° R ip) = G 8 F^(A{R„u^)+u;i(Au,)(R^}, (3.3) 

where R — V, A, denoting the vector and axial-vector resonances respectively. In the 
factorization, the couplings are 

W f = -v« = V2^f RVR , (3.4) 

with r] R is the factorization parameter. f v and Ja are the effective couplings in the 
general strong/electromagnetic lagrangian involving spin-1 resonances (We use the 



notations in Ref. As shown in Refs. [20, 16|, the spin-1 resonance contributions 



to aj's have been obtained. Indeed, the operators in eq. ( p.3|) with R = V do generate 
the structure of the VMD form-factor starting at 0(p 4 ) once the full propagator of the 
vector resonance is taken into account (the axial- vector only contributes a constant). 
The use of the full VMD propagator was also suggested in K — > ttj* |21|] . 

The next leading order magnetic amplitude of K — > Tnr'y in %PT is at 0(p 6 ), 
which contains two parts: local contribution and loop contribution. Although the 
general local 0(p 6 ) couplings in weak effective lagrangian have not been developed 
yet, we can be sure that many unknown parameters have to appear, which will make 
the prediction impossible. One may expect that the local terms could be generated 
through resonance exchange which are reasonably thought as the most relevant ones, 
and the large number of unknown couplings could be reduced significantly. But a 
complete determination of the contributions from resonances including vector, axial- 
vector, scalar, and pseudoscalar still remains very difficult because we have to face 
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some unknown pseudoscalar-resonance weak couplings which cannot be fixed by ex- 
periments. Also, at 0(p 4 ), we know that the couplings a^s have other contributions 
than that from the resonances: for instance, there exists the contribution from the 
WZW anomaly action, giving < a^ n < 1.0 (14|, ^] with a^ n is the unknown param- 
eter. In the Kl case, there is 0(p 6 ) F\ term, which is sensitive to the octet symmetry 
breaking ||, [16|] . 

On the other hand, one can reasonably assume that the photon energy depen- 
dence of the amplitude is dominated by the vector resonance exchange while the 
other resonances including axial-vector, scalar and pseudoscalar only generate the 
constants contributing to the amplitude. Moreover, we have checked the 0(p 6 ) loop 
would lead to very negligible energy dependent contribution in Kl — > 7r + 7r - 7 and 
K + — > 7r + 7T°7. Therefore, a phenomenological description is that, one can express 
the full magnetic amplitude as non-VMD part (which is the constant but with large 
uncertainties involved in it) and VMD part (which is energy dependent and could 
be determined up to one parameter), and use the corresponding experimental de- 
cay rate to determine the former part. Here we use the recent observed values of 
Kl — > 7r + 7r~7 and K + — > 7r + 7r°7: 



Br(K L tt + tt- T , E* > 20MeV) DE = (3.10 ± 0.05) x 10~ 5 [[17 



(3.5) 



and 



Br(K + -> 7r + 7r u 7; 55MeV < T* < 90MeV) DE = (4.7±0.8) x 10" b |23[. (3.6) 



The VMD part magnetic amplitude of Kl — * 7r + 7r 7, corresponding to Figs. 1 
and 2, gives 
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(3.7) 



with 



32V27r 2 /y/i 



v 



(3. 



where hy is the coupling in the general strong/electromagnetic lagrangian involving 
spin-1 resonances fT6f ; the r]v part is 0(p 4 ), and the rest is 0(p 6 ). 
We would like to give some remarks here: 

(1) The VMD form factor like eq. (|1.1| ) was firstly suggested by Lin and Valencia 
[0, and it is phenomenologically successful. However, as already noted by Picciotto 
[p^| , theoretically, there exists some inconsistency in that version because the indirect 
VMD form factor should vanish at 0(p A ) but the 0(p 6 ) indirect vector exchange 
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contribution is important in understanding the magnetic transition of Kl — > 7i + n 7 



24], f|. Here, we have included both direct and indirect VMD contributions to the 
form factor eq. ( |3.7|) with the former one starting at 0(p 4 ) and the latter one starting 
at 0(p 6 ), which properly satisfies all the theoretical constraints. 

2 

(2) The second term in eq. ( |3.7j ), divided by (1 — ^lfz 3 ), is absent in eq. ( |1 . 1| ) . 
However, theoretically, it is generated by Fig. lb (corresponding to 2:3 part) and Fig. 
2b (corresponding to rjv part), we have no reason to exclude it. In fact it construc- 
tively enhances the slope of z 3 , thus affects the rate and the spectrum significantly. 
Note that, as a good approximation, we have used z + = z^ ~ z 3 /2 in deriving this 
term. 

(3) The uj% term in eq. ( |3.3j ) does not contribute to Kl — > 7r + 7r~7, so we do 
not need the factorization relation eq. ( |3.4|) in deriving eq. (|3.7|) . This means the 
present calculation in K L case is independent of the factorization. It is an almost 
model-independent prediction. 

The non- VMD part can be written as 

m l _ eG sm 3 K L , . 

Jw non-VMD — 2lX 2 F ' 

with 

A L = (02 + 2a4) non -vMD + other contributions. (3.10) 



rj v 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



A 0.55 0.36 0.17 -0.01 -0.20 -0.38 -0.57 -0.75 -0.94 -1.13 



A + 1.53 1.34 1.16 0.97 0.79 0.60 0.42 0.23 0.05 -0.13 



Table 1: The quantities A L and A + extracted from the observed branching ratio eqs. 
Q3.5|) and (|3.6|) for the different rjy. 



The values of A L are displayed in the second line of Table 1 in the range of 
0.1 < i]v < 1.0. We find that the spectrum of the photon energy is not very sensitive 
to the value of rjv, as shown in Fig. 3, and they are in agreement with the one 
generated from eq. ( |1 . 1| ) , which gives the best \ 2 fit f° the data in Ref. From 



Fig. 3, the best value of rjv is about 0.5, which is reasonably consistent with the 
preferred rjv — 0.3 obtained in the factorization [|I9| |. The difference is that, some 
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other 0(p 6 ) contributions parametrized using r]yp 7 and ^vpp are considered there 
113, 0. Therefore, it seems that the high order contributions could enhance the 



value of rjy. 



t5 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

E/(GeV) 

Figure 3: E* spectrum of the DE Kl — > tt + it~j (E* > 20MeV). The solid line is from 
eq. (|l.l|) . The rest are generated from eq. (|3.7| ): the dashed line corresponds to rjy = 0.1, 
the dotted line corresponds to rjy = 0.5, and the dot-dashed line corresponds to rjy = 1.0. 



Likewise, in K + — > 7r + 7r°7, the VMD (from Figs. 1 and 2) and non-VMD parts 
magnetic amplitudes are 

eG 8 m 3 K _ 



^VMD 



4tt 2 F 



1 - 2% 



m 



v 



in 



my 



+ 



2 Z 3 



— T~ z + 



Vv + 



+ 



rn 



-Z() 



v 



2m\ 

2~ Z 

my 



and 



with 



^non-VMD 



eG 8 m 3 K 
' 4tt 2 F 



A" 



) 

(3.11) 

(3.12) 
(3.13) 



A + = 2 + 3(2a3 — a2)non-VMD + other contributions. 

We have shown in Table 1 the values of A + for the different rjy, and plotted the 
T* normalized to rrix and W spectrum from the DE magnetic amplitude in Figs. 4 
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and 5, and the W spectrum from the sum of IB and DE amplitude normalized to the 
IB spectrum in Fig. 6. Also,we find that these spectra are not sensitive to the rjy, 



and the last one can be compared with the corresponding experimental result |23| . 



"O 

i 




Figure 4: Spectrum in t*(=T*/m K ) of the DE K + -> 7r + 7r°7 with 55 MeV< T* <90 
MeV. The solid line corresponds to ijy = 0.1. The dashed line corresponds to rjy = 1.0. 



4. Conclusions 



We have presented a phenomenological description of the magnetic amplitudes of 
Kl — > 7r + 7r~7 and K + — > 7r + 7r°7 beyond the leading order in x?T. The VMD 
contribution plays an important role in the analysis. We parameterize the VMD 
part magnetic amplitudes of these two decays in eqs. ( |3.7| ) and fl3.11| ), and the non- 
VMD parts of the amplitudes are estimated by fitting the corresponding observed 
decay rates. Our phenomenological description is consistent with the factorization 
prediction rjy ~ 0.3. 

We summarize the analysis as follows. 

(1) We get the values of A L and A + in the range of 0.1 < rjy < 1.0 (see Table 
1). A L is equals to (ei2 + 2a4) non -vMD plus other higher order contributions. We 
know some other high order contributions, for instance, F\ term is important and 
very sensitive to the octet symmetry breaking [16|]. So here we cannot expect the 



conclusive information on 0,2 + 2a^ from A L . The situation in K + case seems a little 
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O.I 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



W 



Figure 5: Spectrum in W of the DE K + — > 7r + 7r°7. The solid line corresponds to 7?y = 0.1. 
The dashed line corresponds to rjy = 1.0. 

better. After neglecting higher order contributions to A + , we can get —0.71 < (2a 3 — 
a2)non-VMD < —0.1. From Ref. |2U| in the factorization, (2a 3 - a 2 ) axial- vector — 0.3t?a 
with < r]A < 1-0 is the factorization parameter. If we assume the rest contribution 
to (2a 3 — a 2 ) is dominated by the one from WZW anomaly action, we find our 
prediction on (2a 3 — a 2 ) a,n is consistent with the expected < af 1 < 1.0. 

(2) Although from our analysis A L could be positive or negative, the large CP 
asymmtry Bcp |25j in Kl — > 7r + 7r~e + e~ originated from the interference between the 
magnetic and IB amplitude of Kl — > 7r + 7r~7* is predicted to be always positive and 
not sensitive to rjv in the present analysis, which is consistent with the measurement 
p6|. This is not surprising if we carry out Taylor expansion over the form-factor in 



eq. ( |3.7| ) (we assume we can do this expansion), and express the total amplitude as 

eGsm 3 K _ 



M 



-m(l + rz 3 + sz 3 ), 



(4.1) 



2vr 2 F 

we will get, in the range of 0.1 < r/y < 1.0, 

1.47 < m < 1.75, 
2.08 < -r < 2.88, 

2.50 < s < 3.93, (4.2) 

which are comparable with the recent KTeV measurement r = —2.93 ± 0.41 ± 0.34, 
s = 3.31 ± 1.15 ± 0.96 [17], and |m| = 1.53 ± 0.25 in Ref. from only linear slope 
fit. 
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2 - 
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0.5 - 



0.1 




Figure 6: W spectrum normalized to the IB spectrum of K + — ► 7r + 7r°7. The dashed line 
corresponds to r\y = 0.1. The dotted line corresponds to r/y = 1.0. The data points are 
fromRef. Wc 



(3) So far, there is no experimental evidence for the energy dependence of the 
magnetic amplitude in K + — ► 7r + 7r°7. But the VMD form-factor obviously indicates 
this energy dependence. By Taylor expansion of eq. ( ft.ll ), we can get the corre- 
sponding linear slopes of z + and z$: 0.77 < — r + < 1.61, 0.72 < — r < 1.17 in the 
range of 0.1 < rjy < 1.0. These values are not small if this kind of Taylor expansion is 
valid here. Unfortunately, it is not very easy to measure these quantities experimen- 
tally. On the other hand, we find that, Zq is related to T* through a linear relation 
eq. ( p.6| ). Therefore, it is expected that a high-precision experimental analysis of 
the T*/rrtK distribution from DE contribution may be able to measure this energy 
dependence of the amplitude. 
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